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system  of  interest  is  the  lower  order  observer  or  lower  order  estimator. 

The  determination  of  changes  in  the  performance  of  the  lower  order 
observer  to  changes  in  parameters  is  of  great  importance  in  engineering 
analysis  and  design.  Given  the  fact  that  the  output  of  the  lower  order 
observers  are  physically  not  accessible,  a study  of  the  changes  of  these 
unaccessible  states  to  changes  in  system  parameters  are  helpful  in 
' determining  the  properties  of  a system.  Besides  this,  observers  are 

■ used  for  adaptive  schemes  in  which  the  sensitivity  functions  are  needed. 

I 

The  significance  of  using  lower  order  observers  instead  of  a full  order 
observer,  when  a few  states  are  needed  to  be  estimated,  is  of  an  economical 
< nature.  Sometimes  it  is  needed  that  the  estimated  state  be  highly  non- 
' sensitive  to  parameter  variations,  the  economic  reasons  are  then  super- 

’ seded  by  sensitivity.  Therefore,  a comparison  in  sensitivity  is  needed 

between  the  lower  order,  and  full  order,  observers. 
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1.  INTRODUCTION 

The  problem  of  sensitivity  of  dynamic  systems  to  changes  in 
parameters  has  acquired  great  importance  in  modern  technology.  Sensitivity 
is  important  because  in  physical  realization  of  systems  of  automatic  control, 
changes  in  the  parameters  always  are  encountered  as  a result  of  aging  of  the 
elements,  the  effects  of  the  external  environment,  interactions  with  other 
systems  etc.  Furthermore,  contemporary  automatic  control  systems  often  are 
realized  as  systems  of  variable  structure  where  the  possibility  of  varying 
the  parameters  of  the  system  is  specially  introduced  to  obtain  adaptation 
properties  as,  for  example,  in  systems  of  optimal  control  or  in  self-learning 
systems . 

In  all  these  cases  the  question  arises  as  to  how  the  behavior  of 
the  system  changes  when  changes  occur  in  certain  of  its  parameters.  One 
such  system  of  interest  is  the  lower  order  observer  or  lower  order  estimator. 

The  determination  of  changes  in  the  performance  of  the  lower  order 
observer  to  changes  in  parameters  is  of  great  importance  in  engineering 
analysis  and  design.  Given  the  fact  that  the  output  of  the  lower  order 
observers  are  physically  not  accessible,  a study  of  the  changes  of  these 
unaccessible  states  to  changes  in  system  parameters  are  helpful  in  determin- 
ing the  properties  of  a system.  Besides  this,  observers  are  used  for  adaptive 
schemes  in  which  the  sensitivity  functions  are  needed. 

The  significance  of  using  lower  order  observers  instead  of  a full 
order  observer,  when  a few  states  are  needed  to  be  estimated,  is  of  an 


economical  nature.  Sometimes  it  is  needed  that  the  estimated  state  be 
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highly  nonsensitive  to  parameter  variations,  the  economic  reasons  are  then 
superaeded  by  sensitivity.  Therefore  a comparison  in  sensitivity  is 
needed  between  the  lower  order,  and  full  order,  observers. 

In  order  to  study  this  behavior  the  sensitivity  functions  or  the 
parameter  influence  coefficients,  of  the  estimated  states  of  the  system  with 
respect  to  system  parameters,  will  have  to  be  generated.  If  the  system  is 
directly  differentiated  with  respect  to  the  parameters  as  in  [9],  resulting  in 
differential  equations  of  very  high  dimensions.  To  avoid  this  method,  the 
method  of  generating  sensitivity  functions  as  proposed  by  Wilkie  and 
Perkins  in  [ll  2]  is  adapted  to  reduce  the  model  of  the  system. 

Chapter  two  deals  entirely  with  the  development  of  the  sensitivity 
model.  As  an  example,  the  Luenberger  Observer  has  been  used.  The  low  ord- 
er sensitivity  model  is  determined  and  certain  schemes  have  been  set  up 
to  help  design  this  model.  Chapter  three  deals  with  the  experimentation  of 
some  practical  examples.  The  examples  as,  for  example  a d-c  motor  and  a 
thermometer  are  used  in  order  to  give  the  reader  a physical  feeling  for  the 
analysis  of  these  systems.  Chapter  four  deals  with  the  conclusions  drawn 
from  the  experimental  results  of  Chapter  three. 

Finally,  this  analysis  will  shed  some  light  on  the  properties  of 
the  lower  order  observers.  In  no  way  are  the  conclusions  in  this  thesis 
applicable  to  each  and  every  lower  order  observer,  but  it  is  assumed  that 
the  conclusions  drawn  in  this  thesis  can  be  attributed  to  lower  order 


observers  in  general. 
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2.  SENSITIVITY  tIODEL  OF  THE  LOWER  ORDER  OBSERVER 

2.1.  Introduction 

In  this  chapter,  the  sensitivity  model  for  the  lower  order 
observer  will  be  developed.  It  is  of  considerable  practical  and  theoretical 
interest  to  obtain  the  trajectory  sensitivity  functions.  In  the  case  of 
this  study,  the  functions  will  be  used  to  study  the  advantages  or  disad- 
vantages of  the  lower  order  observer,  which  is  carried  out  in  Chapter  3. 

From  a practical  point  of  view,  low-order  realizations  are  required  for 
efficient  use  of  equipment  in  simulations  and  for  reasons  of  accuracy. 

In  this  chapter  we  shall  show  that  all  the  sensitivity  functions 
can  be  generated  from  a n-dimensional  sensitivity  model,  or  under  certain 
conditions,  from  a (n-r)-dimensional  model. 

An  overview  of  sensitivity  is  available  in  [3,4],  The  sensitivity 
model  developed  in  this  chapter  is  a direct  application  of  the  model  pro- 
posed by  Wilkie  and  Perkins  in  [ll,12], 

2.2,  Lower  Order  Observers 

The  state  variables  of  a system  are  sometimes  needed  for  various 
purposes,  the  common  one  being  state  feedback.  When  designing  a feedback, 
it  is  desired  that  all  the  states  of  the  system  are  available.  This 
desired  situation  often  does  not  hold  in  practice,  however,  either  because 
the  state  variables  are  not  accessible  for  direct  measurement  or  because  the 
number  of  measuring  devices  is  limited.  Therefore,  in  order  to  stabilize, 
decouple  and  to  optimize  a system  that  has  a limited  number  of  states 
available,  a reasonable  substitute  for  the  state  vector  has  to  be  found. 

The  dynamical  system  that  constructs  an  approximation  of  the  state  vector 
Is  called  a state  estimator  or  an  observer. 
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There  are  many  different  kinds  of  observers,  A textbook  summary 

of  this  subject  is  given  in  [2].  Generally  there  are  two  types  of  observers. 

The  n-dlmensional  or  full  order  observer  and  the  lower  order  observer.  The 
n-d imens ional  full  order  observer  estimates  all  the  states  of  the  system. 
Sometimes  this  is  unnecessary,  as  some  of  the  states  may  be  obtained  exactly 
from  measured  outputs.  In  such  cases  a lower  order  observer  is  used  to 
estimate  only  those  states  that  are  required. 

The  lower  order  observer  can  be  designed  in  many  ways,  A text- 
book treatment  of  observers  may  be  found  in  [2],  In  this  study  we  shall 

consider  the  Luenberger  Observer  [8, 13]  . We  shall  use  the  circumflex  over  a 

variable  to  denote  an  estimate  of  the  variable.  For  example,  x^  is  an 
estimate  of  x^;  x^  is  an  estimate  of 

Consider  the  time-invariant,  n-dimensional  system  described  by 
the  state  equations 

x=Ax+Bu  (2.2.1) 

y =•  Cx  (2.2.2) 

The  output  y is  r -dimens ional. 

The  matrices  A,  B and  C are  assumed  to  be  known.  The  problem  is 
to  generate  all  state  variables  from  the  available  input  u and  output  y. 

The  equation  of  the  Luenberger  Observer  as  given  in  [5,7]  are 

Z - PZ  + 0y  + Ru  (2,2.3) 

x^  = Z + Ky  (2.2.4) 
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where 

(n-r)  X 1 state  vector 
Z is  the  (n-r)  X 1 state  vector 
P is  an  (n-r)  X (n-r)  real  constant  matrix 
0 is  an  (n-r)  x r real  constant  matrix 

R is  an  (n-r)  x n real  constant  matrix 

K is  an  (n-r)  x r real  constant  matrix 

Assume  C has  rank  r.  There  is  no  loss  in  generality  in  this 

assumption.  If  rank  C < r,  some  outputs  are  linearly  dependent  on  others 
and  can  be  eliminated.  Then  there  exists  an  (n-r)xn  matrix  such  that 
the  transformation  N in 


is  non-singular. 


X 

a 

A A 

= Nx  = 

‘S 

C 

- . 

x 


Then 


NAN 


-1 


11 


‘12 


NB 


(2.2,5) 


(2.2.6) 


Therefore 


X 

a 


A X 
11  a 


+ A^2^S  ^l'" 


(2.2.7) 


jc  = A x^  + A„3c  + B,u 


(2.2.8) 
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The  appropriate  choices  of  matrices  for  (2,2.3)  and  (2.2.4)  can 
then  be  shovm  to  be 


P = A^l  - KA2^  (2.2.9) 

0 = (A^^  - KA2^)K  + (A^2  • <^22^  (2,2.10) 

R » - KB2  (2.2.11) 

The  matrix  K is  selected  so  that  the  eigenvalues  of  Aj^j^  - KA2j^ 
are  the  desired  eigenvalues,  and  make  the  error,  decay  to  zero.  Figure  2,1 
is  the  block  diagram  of  the  Luenberger  Observer. 

2.3.  Sensitivity 

The  system  given  by  (2.2.1)  can  be  dependent  on  many  parameters. 

In  general  let  the  parameter  vector  effecting  system  (2.2.1)  be  given  by  g, 
where  vector  g is  N-dimenslonal . 

A variation  in  the  parameter  g will  produce  a variation  in  the 
solution,  which,  to  first  order,  is  given  by 


6x^(t,g) 


N 


dg 


j' 


The  partial  derivatives 


j = 1,2,... ,N  (2.3.1) 


^ A 


I - 1,2 n-r 

1 - 1,2,. ..,N 


(2.3.2) 


are  called  the  estimated  state  trajectory  sensitivity  functions. 
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These  sensitivity  functions  are  of  prime  importance  in  determining 
the  effect  of  parameter  variations  on  the  effectiveness  of  the  observer 
in  estimating  the  state. 

The  linear  time-invariant,  n-dimensional  system  described  by 
(2.2.1)  is  now  written  as 

x(g)  = A(g)x(g)  + B(g)u  (2.3,3) 

y(g)  = C x(g)  (2.3,4) 

where  g is  a N-dimensional  parameter  vector. 

Considering  the  lower  order  observer  given  by  (2.2,3)  and  (2.2.4), 
it  must  be  noted  that  the  matrices  P,  9,  R and  K are  constants  and  are 
desired  at  a nominal  value  of  the  parameter  g,  given  as  g^. 

Thus  the  functional  dependence  of  (2.2.3)  and  (2.2,4)on  g is 

Z(g)  = PZ(g)  + 0y(g)  + Ru  (2.3.5) 

x(g)  = Z(g)  + Ky(g)  (2.3.6) 

2,4.  Sensitivity  Model  for  the  Lower  Order  Observer  in  Open-loop 

In  the  introduction  to  this  chapter,  the  significance  of  reduced 
models  was  explained.  The  sensitivity  functions  for  linear  tlme-1 nvariant 
systems  may  be  obtained,  in  principle,  by  solving  the  linear  equations  with 
literal  parameters  and  differentiating  the  result.  This  procedure  is 
extremely  tedious  and  un illuminating.  A more  useful  approach,  as  shown  by 
Wilkie  and  Perkins  in  [11,12],  is  to  obtain  differential  equations  for 
the  trajectory  sensitivity  functions.  Furthermore,  the  structure  of  the 


resulting  simulation  itself  often  provides  helpful  Insight  into  the  system 
behavior  with  respect  to  the  parameters.  The  Wilkle-Perkins  low  order 
sensitivity  model  will  now  be  developed  for  the  observer  given  by  (2.2.3) 
and  (2,2.4),  It  must  be  noted  that  the  sensitivity  functions  needed  are 
the  estimated  state  sensitivity  functions  only  and  not  the  sensitivity 
functions  of  the  available  states.  We  shall  now  consider  the  single  input 
case  of  system  (2.2.1), 

Now  consider  (2.2.7).  It  is  a (n-r )-dimensional  equation  and  is 
the  differential  equation  of  all  of  the  states  that  have  to  be  estimated. 
Then,  the  sensitivity  model  for  only  the  lower  order  observer  will  be  a 
(n-r )-dimensional  model.  It  is  seen  that  the  input  to  this  model  are  the 
sensitivity  functions  of  the  output  of  system  (2.2.1),  therefore,  we  will 
be  required  to  generate  the  output  sensitivity  functions  of  system  (2.2.1). 

In  order  to  accomplish  this,  we  substitute  for  in  (2,2,5) 


and  get 


X 


(2.4.1) 


The  vector  is  similar  to  the  state  vector  x and  therefore  we 


-ir.-.. 


I 
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where  b is  a n x 1 matrix.  Note  that  system  (2. A. 2),  (2.4.3)  has  matrix 
A which  is  similar  to  system  (2.2.1)  and  therefore  (2.4.2),  (2.4.3)  is 
considered  to  be  an  open-loop  system.  We  shall  now  proceed  to  develope 
an  open-loop  sensitivity  model  for  (2.4.2),  (2.4.3). 

In  [llj  it  has  been  shown  that  for  the  controllable  system 
(2.4.2)  and  (2.4.3)  a non-singular  transformation 


X = Tq 


(2.4.4) 


always  exists,  such  that 


q = Aq  + bu 


(2.4.5) 


where 


A = t"^AT 


0 

1 

0 . . 

. o" 

0 

0 

0 

1 . . 

. 0 

0 

• 

. 

• • • 

. . 

1 

H 

II 

1-0 

• 

0 

0 

0 . . 

. 1 

• 

^1 

~°‘2 

-Qf^  . . 

.-a 

n 

1 

L «i 

(2.4.6) 


The  matrix  T can  easily  be  found  as  shown  in  [3] . The  elements 
Qf^  ,©2  f • • • i*'  matrix  A are  the  coefficients  in  the  characteristic  equation. 

It  is  evident  that  A,  b and  T are  a function  of  the  parameter 
vector  g.  The  sensitivity  functions  for  these  parameters  will  now  be 
found. 

It  can  easily  be  shown  that  (2.4.5)  can  be  written  as 


qCt.a)  - A(o')q(t,a)  + bu(t) 


(2.4.7) 


- 


! 


W 


P 

I 
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Taking  the  partial  deravative  of  both  sides  of  (2.4,7)  with  respect 
to  O'  and  interchanging  the  order  of  the  differential  equation  we  get 


6a 


6a  ^ 

1 = 1 ,2 , . . . ,n 


Let 


6a 


, i=l,2,.,.,n 


(2.4.8) 


(2.4.9) 


such  that  the  initial  conditions  are 


hj  = h.(t,a)|^^Q  . 


(2.4.10) 


Assuming  appropriate  continuity  conditions  we  Interchange  the 
order  of  the  partial  with  respect  to  A and  substitute  t=0. 

Then, 


^0  ” 6a, 


6a, 


(2.4.11) 


From  (2.4.  8 ) and  (2.4.  9 ) we  get 


hj^(t.a)  = A(a)h^(t,a)  + q(t,a) 


(2.4.12) 


From  the  differential  equations  given  by  (2.4.12)  the  sensitivity 
functions  for  the  estimated  states  can  be  easily  found.  The  subject  is 
studied  thoroughly  in  Section  2.5. 

There  are  two  Important  properties  of  the  sensitivity  model.  As 
the  model  is  similar  to  the  one  described  by  Wilkie  and  Perkins  in  [11,12] 
the  proof  of  these  properties  will  be  omitted,  but  the  properties  will  be 
stated  for  convenience. 
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Total  Symmetry  Property 

The  sensitivity  functions  defined  by  (2.4.9)  is  given  as 


i 


and  a sensitivity  matrix  as 


[h  .]  = [h  h h„  . . , h ] 
1 9 j 1.  £.  i.  n 


(2.4.13) 


Th 


property 


e sensitivity  matrix  [h.  .]  has  the  following  total  symmetry 

^ » J 


1 ^ “ 1,2,. ...n  (2.4.14) 

i-i.J+i  j = 1.2,...,n-l 


Then  all  the  elements  along  the  antidiagonals  of  [h  ] are  equal 

i > J 

as  shown  in  (2.4.15). 


1.1 

^.2 

1,3  l,n 

1,2 

- 

|'2.n 

1.3 

h 1 

n-1  ,n 

l,n 

‘^2,n 

• • • h • h 

n-l,n  n,n 

(2.4.15) 


Complete  Simultaneity  Property 

All  the  sensitivity  functions  , i,  j “ l,2,...,n,  can  be 
obtained  as  linear  combinations  of  the  states  of  the  sensitivity  model 
(2.4.12)  and  the  system  states  q. 


r* 

i 


I 
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Thus  far  we  have  developed  a n-dimensional  sensitivity  model. 

The  important  point  of  this  model  is  the  fact  that  all  the  sensitivity 

functions  can  be  generated  by  this  model.  If  for  example  we  have  j parameters 

in  the  system,  then  all  the  (n-r)  X j sensitivity  functions  can  be  generated 

by  the  n-dimensional  model.  It  is  important  to  state  at  this  point  that, 

6A 

by  carefully  choosing  the  input  to  the  sensitivity  model  can  be  an 

available  state  of  system  (2.2.1). 


2.5.  Sensitivity  Functions 

In  Section  2.4  a sensitivity  model  was  developed  for  the  observer, 
which  is  given  by  (2.4.12).  We  will  now  show  how  the  estimated  state 
sensitivity  functions  will  be  generated  from  one  model  of  n-dimension. 

From  the  transformations  (2.2.5)  and  (2.4.4)  we  have 


X = N,Tq 
a 1 


(2.5.1) 


Note  that  T is  a function  of  the  parameter  vector  g. 

Taking  the  partial  of  both  sides  with  respect  to  a we  get 


^(t,a) 

6a  ^ 


From  (2.3.2)  we  have 


6[N  T(g)] 

[N^T]h^(t,a)  + ^ q(t,a) 


Ji 

i 


1,2,. 

1,2,. 


,n-r 

.n 


(2.5.2) 


6x 


a£  A - 
6gj  “ 


jI  - 1,2, .. . ,n-r 
j - 1.2 N 


(2.5.3) 


f 


60  . 6 g . 

1 


. ,n-r 
.,N 


(2.5.4) 


Substituting  (2.5.2)  in  (2.5.4)  we  have 

6[n  tJ  6a. 

{[NiTlh  (t.a)  + q(t.a)].^  . C 

1 Dj  .J 

i * 1,2, ... ,n 
I « 1,2, .. . ,n-r 
j = 1,2,...,N 


(2.5.5) 


As  C is  independent  of  g,  then  Nj^  is  also  independent  of  g.  We 
6a  6a 

can  interchange  h 7 — with  7 — h.  Therefore 
og  6g 


—O 


(2.5.6) 


i » 1,2, ...  ,n 
I = 1,2,... ,n-r 
j = 1,2,... ,N 

The  equation  (2.5.6)  gives  all  the  sensitivity  functions,  (n-r)  x N. 
Thus  we  see  that  a n-dlmensional  nuidel  gives  all  the  (n-r)  X N estimated 
state  sensitivity  functions.  It  is  very  easy  to  show  how  these  functions  are 
generated.  The  states  h (t,a)  are  derived  from  the  sensitivity  model  (2.4.12), 

1-  e _ 


q(t,a)  is  generated  by  transformation  (2.4.4).  The  J easily 

calculated  as  the  functions  a^  and  are  known  and  can  be  easily  calculated 
on  a computer.  For  example  If  gj^  * a^  then  (2.5.6)  will  be  given  as 


> - c, 


. ,n-r 
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In  tho  same  way  all  the  sensitivity  functions  can  be  found.  It  is  very 

important  to  note  here  that  in  the  case  of  single  input,  the  maximum  number 

2 

of  parameters  are  n +n. 

It  has  been  shown  that  the  matrix  T can  be  calculated  by  the 
relations  1 , 14J 

T = [t,  t„  . . . t ] 12. 5, g 

1 Z n 


whore 


t = b 
n 

t , = At  + O'  t 
n-1  n n n 

t « = At  - + cy  , t 
n-2  n-1  n-1  n 


tl  “'z'^n 


Differentiating  (2,5.9)  with  respect  to  we  get 


n ^ 6b 
Sg  j 


u u - 

n-1 

- . . 

n 

figj 

f t 1 

figj  n 

fig. 

'Vi 

figj 

figj  '"n-l 

+ A 

figj 

ftgj 

t-  + A 
figj  2 

64 

figj 

60^2 

n ^ fib 
fig.  n fig^ 


n-1  fig. 


b + O' 


2 fig, 


lb 


Tims,  if  the  functional  dependence  of  the  system  matrix  A is 

^ A n ( 8 ) 

assumed  known  (so  that  j is  known)  and  also  |j i is  known.  The  - — ^ ' 

Its, I ->gj 

can  be  easily  calculated.  Experience  has  shown  that  when  using  computers 
the  time  needed  for  the  algebraic  calculations  associated  with  the  trans- 
formation to  canonic  form  is  much  less  than  the  time  required  for  the 
additional  integrations  needed  in  direct  implementation  of  the  model. 

Figure  2.2  is  the  block  diagram  of  the  sensitivity  model,  as  proposed  in  Sec- 
tion 2.5.  A few  examples  have  been  worked  out  in  Chapter  3. 


2.6.  Sensitivity  Model  for  the  Lower  Order  Observer  in  a Feedback  Loop 

It  has  been  mentioned  before  that  the  lower  order  observer  is  used 

frequently  for  state  feedback.  It  would  be  appropriate,  then  to  develop  a 

model  for  the  closed-loop  sensitivity  functions.  Using  the  same  technique 

as  developed  in  Section  2.5,  we  will  now  show  liow  to  develop  this  model. 

Consider  (2.4.7).  Let  the  feedback  be  given  by 

u = -u-(-u.  (2.6.1) 

a b 

where 

u^  = Cq  (2.6.2) 

Let  matrix  <t  = [(J,  (S-  . . . (}  ] = FT,  where  F is  a constant  matrix  and  T 

L z 3 n 

is  the  transformation  matrix  from  (2.4.4),  dependent  on  the  parameter 
vectors  a from  A and  6 from  matrix  b.  The  input  u^  is  then,  a function 
of  the  parameter  vectors  a and  3.  Let  the  input  u^  be  a function  of  time. 
Note  that  the  coefficients  of  the  parameter  vector  B are  obtained  from  the 
matrix  b.  A textbook  treatment  of  state  feedback  pole  placement  may  be  found 
in  TlO]. 


Combining  (2,4.7)  and  (2.6.1)  we  get 


q(t,a,e)  = (A(a)  - b(j5(Q',3)jq(t,Q',3)  + buj^(t) 
Equation  (2.6.3)  can  then  be  written  as 


(2 


where 


and 


q(t,uj)  - W(uD)q(t,uu)  + buj^(t) 


W = 


0 1 0 

0 0 1 

• • • 

0 0 0 


-0J2 


ou^  = 

(^2  = 0^2  ®2 
0)3  - 0^3  + 03 


(2 


(2 


(2. 


The  parameters  u)j^,uu2,. 


Let 


(JU  “O'  +0 

n n n 


,U)^  are  the  coefficients  of  the  parameter  vector 


M(t.u))  A 

6u). 


(2. 


.6.3) 


.6,4) 


.6.5) 


6.6) 


U). 

6.7) 
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such  that  the  initial  conditions  are 


s'  = S (t  U))l  = §illLr!£il 
^0  i*-  ’ '’'t=0  6UU.  't=0 


(2.6.8) 


also 


Mr 


i 6 , , . I , ^0 

1 1 


(2.6.9) 

Using  the  same  techniques  as  used  in  Section  2 . 5 we  can  easily 


show  tliat 


S.(t,x)  = W((r)S^(t,UJ)  +^|^q(t,o)).  (2.6.10) 

i 

Equations  (2.6.10)  and  (2.4.12)  are  alike  excepting  that,  (2.6.10) 
is  dependent  on  a parameter  vector  uj  and  (2.4.12)  on  a parameter  vector  a. 
Similarly,  it  can  be  shown  that  (2.6.10)  also  has  the  two  properties  of 
(2.4.12)  and  therefore  we  can  generate  all  the  sensitivity  functions  from 
this  n-d imens iona 1 model. 

I'he  equations  of  the  sensitivity  functions  are  then  given  as 


(2.6.11) 


i = 1 , 2 , . . . , n 
Z = 1,2,..., (n-r ) 
i = 1 , 2 , . . . , r 

These  are  the  equations  of  the  sensitivity  functions.  It  must 
be  noted  that  the  closed-loop  model  and  the  open-loop  model  are  similar, 
excepting  that,  whereas  the  open- loop  model  has  parameter  vector  a,  the 
closed-loop  model  has  parameter  vector  u;.  Using  the  same  discussion  as 
before  we  can  easily  compute  I | ’ computed  as  shown  before. 

This  concludes  our  discussion  of  the  sensitivity  model.  In  the 
next  chapter  we  sliall  generate  the  sensitivity  functions  of  a few  examples 
and  analyze  them. 


20 


3.  ANALYSIS  OF  THE  SENSITIVITY  OF  LOWER  ORDER  OBSERI'ERS 
3.1.  Introduction 

When  designing  a lower  order  observer  it  is  very  important  to  know 
how  good  the  estimation  will  be.  One  of  the  fundamental  requirements  is 
to  reduce  ttic  error  between  the  output  of  the  system  and  the  estimated  states, 
finally  decaying  to  zero  in  finite  time.  In  Section  2.2  we  have  seen  how  this 
is  achieved  by  selecting  the  matrix  K,  in  other  words  by  selecting  the 
eigenvalues  of  the  lower  ordef  observer.  -This  approach  to  the  problem  is 
fine,  assuming  that  notliing  else  effects  the  estimated  states.  In  practice, 
however  this  is  not  true.  One  of  the  reasons  is  that  the  parameters  of  a 
system  are  prone  to  variations,  thus  affecting  tlie  estimated  states.  The 
error  then,  is  no  longer  only  dependent  on  the  eigenvalues  of  the  lower 
order  observer  but  also  on  the  parameters  of  the  system.  If  the  error  for  a 
particular  estimated  state  is  required  to  remain  within  fixed  limits  and  to 
decay  to  zero  in  finite  time  then  there  may  be  a set  of  eigenvalues  of  the 
lower  order  observer  that  will  accomplish  this.  But  there  may  be  one  single 
parameter  in  the  system  that  could  drive  the  error  beyond  the  limits,  thus 
eliminating  the  possibility  of  using  that  particular  lower  order  observer. 

It  lias  been  mentioned  before  that  lower  order  observers  are  used 
extensively  for  the  purpose  of  state  feedback.  Many  investigators  have 
mentioned,  that  the  sensitivity  of  a system  improves  when  using  feedback. 

I t is  therefore  intriguing  to  study  the  sensitivity  of  lower  order  observers 
in  a feedback  loop.  Another  aspect  that  is  worthy  to  study  is  that,  given 
a set  of  allowable  eigenvalues  of  the  lower  order  observer,  which  eigenvalues 
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uill  c .lusi  Lliu  lower  order  observer  to  be  least  dependent  on  the  parameters, 
that  is,  eflects  iif  the  eigenvalues  oC  a lower  order  observer  on  its 
sensitivity.  I'his  is  a very  important  element  in  optimizing  the  best 
eigenvalues  in  order  to  design  the  best  observer. 

It  was  mentioned  earlier  that  one  of  the  reasons  that  lower  order 
obsi  rvers  weie  preferred  to  full  order  observers  was  economics.  Sometimes 
this  could  i>rove  untrue,  when  sensitivity  supersedes  economics  and  a highly 
i nscii.s  i t i ve  observer  has  to  be  designed.  The  question  then  arises,  of  how 
do  lover  order  observers  and  full  order  observers  compare  in  sensitivity. 

This  is  another  point  that  needs  investigation. 

Thus,  we  see  that  a study  of  the  sensitivity  functions  of  the 
estimated  states  will  give  us  an  insight  into  some  of  the  aspects  of 
designing  a bitter  lower  order  observer  or  replacing  i t by  a full  order 
ob.se.rver.  In  this  chapter  we  shall  investigate  the  above  mentioned  points. 
S>)me  practic.il  e..amples  will  be  used  to  illustrate  the  subject  matter.  It 
mu.'.t  be  nuieil  here  that  the  analysis  will  cover  only  the  topics  mentioned, 
eviii  th.ough  there  may  be  many  more  related  aspects  that  need  to  be 
1 nvest i gated . 

Examples  I,  2,  3,  and  4 generate  the  open-loop  sensitivity  functions 
and  thus  will  be  used  for  open-loop  sensitivity  analysis.  Examples  5,  6,  and 
7 are  used  for  the  closed-loop  analysis  of  the  lower  order  observer  sensi- 


tivity functions. 
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3.2,  Example  No.  1 

The  equations  of  an  armature  coupled  d-c  motor  is  given  by  the 

equations 


where 


L , LB+RJ«  RB  + K, 

( jn  js^  + b T 3 0 

Kt  ' m 


9 is  tlie  angle  of  displacement  of  the  motor 
m 


e input  voltage  to  the  armature 

2i 


L inductance  of  the  armature 


m 

R resistance  of  the  armature 


m 

J moment  of  inertia  of  motor  plus  load 
B damping  coefficient 

and  proportionality  constants. 


(3.2.1) 


The  nominal  value  of  these  parameters  are  given  by 

L°  “ 1 henry 
m 

R°  “ 2 ohms 
m 

- 1.125  Kg-m^ 

B°  “ 1.125  newton/ (meter /second) 

K°  “ 2.25  newton-m/amp  or  volt-sec/radian 
Kf  ■ 1 volts-meter-sec/weber-rad. 


(3.2.2) 


Equation  (3.2.1)  can  be  put  into  the  form 
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0 

1 

0 

' 0 

X = 

0 

0 

1 

X + 

0 

0 

-a 

0 

■“3 

y = x, 


where 


o ^ W _ 2 X 1.125  X 1 X 2.25 

2 L J “ 1.125  X 1 

m 


(3.2.3) 


(3.2.4) 


T R T 

o ^ m _ 1.125  x 1 + 1.125  X 2 

3 L J ~ 1.125  X 1 

m 


gO  _ ^ 2.25 

3 “ ^0^0  “ 1.125  X 1 


(3.2.5) 


(3.2.6) 


The  eigenvalues  of  the  system  are  0,  > "2 

2 2 2 

There  are  six  parameters  in  this  system. 


- B (3.2.7) 

S “ \ 

'’e  ■ 


- y/i . 

2 

They  are 


The  reason  why  these  parameters  may  change  a little  can  be  explained 
when  the  physical  aspects  of  this  type  of  motor  Is  studied  In  detail.  As  for 
now  we  assume  that  these  parameters  are  subject  to  changes. 
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It  is  required  that  states  and  be  estimated.  Therefore 
designing  a l.uenberger  observer  for  system  (2.2.3)  we  have 


z + 


2 1 
- kf  + k 

0 

1 2 

y + 

_»3_ 

(3.2.8) 


The  coefficients  and  are  selected  such  that  the  eigenvalues 
of  (3.2.8)  are  more  negative  than  the  eigenvalues  of  system  (3.2.3). 

The  output  of  the  observer  is 


><2  “ Zf  + k^y  , ^3  “ ^2  + k^y 


(3.2.9) 


With  the  help  of  the  theory  developed  in  chapter  2,  the  open-loop 
sensitivity  model  is  given  by 


~0 

1 

0 

-1 

h “ 

0 

0 

1 

h - 8° 

^^2 

(3.2. 

10) 

0 

o 

-^2 

- 

a° 

3_ 

The  sensitivity  functions  can  now 

be 

obtained 

from  the  model 

(3.2.10).  The 

transformation  matrices  are 

given 

as 

"1 

0 

0“ 

”o  1 ol 

T - 

^3 

0 

1 

0 

> 

^ * 

0 0 1 

(3.2. 

11) 

_0 

0 

1_ 

- J 

As  there  are  two  estimated  states  and  six  parameters,  there  will 
be  twelve  sensitivity  functions  of  the  estimated  states.  They  are  as  follows: 
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. ^ o o 

!!i.  Vi-u  i_- 

5l  ■ 3 o ''2,3  L "2 

m L J ra 

m 


"^^2.3  • "2 


(3,2,12) 


I 

I 


c'*  00 

S_  JL^ 

■L  2 ^ ^3,3  L ^^3 

m L°  jO  ’ 
m 


'■'■3.3  ■ “s 


6x.  k; 


5r 


o2  o 
L J 
m 


*^2,3  “ ^ *^2,3 


(3.2.13) 


(3.2. K) 


i i 


I 


«3 


h.  = 2 h„ 


f>R  ,2,0  “3,3  3,3 

m Ij  j 


(3.2.15) 


‘’‘2  V ^ 1 - 

‘J  ■ o3 

J i.4 


l,78h^  „ - 0,89x^ 
2,3  2 


(3.2,16) 


- ^ U 

'■-3  V . 1 - 


fij 


- h_  . - — x„  = - 1.78h-  „ - 0.89X.J 

3 3,3,o3  3,3  3 

J°  L° 
m 


(3.2.17) 


£ ^ O 

2 T. 

&B 


L°J° 

m 


2 ^2,3  “ ^'^®'^2,3 


(3,2.18) 


c ^ L 

6x„  K„ 


&B 


L°  J° 
m 


h-  - - 1.78h 
2 3,3  3,3 


(3.2.19) 


^^2  « 1 
T ,0  ,0  ' 


L J 
m 


(3.2.20) 
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o o 

' 

2 2 ^2,3  ^ o 3 

L°  J°  "t 

m 


1.78h2  3 + 0.4X3 


1,3 


4 h 


1,3 


(3.2.21) 


(3.2.22) 


4 h 


2,3 


(3.2.23) 


Thus  we  see  that,  with  the  knowledge  of  only  three  sensitivity  functions 
h,  h -,  and  h , we  have  generated  all  the  twelve  sensitivity 
functions,  (3.2.12)  to  (3.2.23).  Figure  3.1  is  the  diagram  for  generating 
the  estimated  state  sensitivity  functions  for  Example  1.  The  plots  of  the 
sensitivity  functions  (3.2.12)  to  (3.2.23)  are  illustrated  in  Figures  3.2 
to  3.5. 

Comments 

While  developing  the  sensitivity  model,  it  is  noticed  that  the 

partitioned  form  of  (3.2.3)  has  A,,  I’e  a zero  matrix.  This  eliminates  the 

6x1 

need  for  the  output  sensitivity  function  and  thus  reduces  the  sensitivity 

model  to  a (n-r)  dimensional  model.  Thus  it  can  be  concluded  that  if 
equation  (2.2.7)  has  rank  (n-r),  then  the  sensitivity  model  for  the  observer 
will  be  a (n-r)  dimensional  model.  If  the  rank  is  higher,  then  we  must 
use  the  n-dlmensional  model  as  derived  in  Section  2.5. 

It  must  be  noted  that  the  model  generates  the  open-loop,  estimated 
state  sensitivity  functions. 


I 

I 
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gure  3.1  Diagram  for  generating  tli«  estimated  state  sensitivity  functions  for  Examp 


I 

I 

I 

I 


0 2 


Time,  in  sec 

Figure  3.2  Estimated  state  >(j_and  t 
functions  for  Example  1. 


I 
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Time,  in  sec 

Flgvre  3.3  Estimated  state  1?j.and  the  sensitivity 
functions  for  Example  1. 


L’ 


It  is  noticed  from  Figures  3.2  to  3.5  that  all  the  sensitivity 
functions  have  greater  magnitude  than  the  actual  estimated  states  and 
An  example  to  what  can  happen  is  given  as  follows. 

For  n ten  percent  increase  in  parameters  the  deviations  in  are 
tabulated  below.  These  deviations  are  the  steady  state  deviations  and 


will  effect  the  steady  state  response  of  . 


Sx^  will  be  1.9 
” ” 1.0125 

" " -1.8 

" " 0.65 

" " 0.675 

" " -0.84375 

6x2  will  be  1.583 


for  a ten  percent  increase  in  L 

n 

M „ M j 

M It  n 


as  a combined  effect. 


Therefore  limit  X2  tnay  always  be  in  steady  state.  This  could  be  a very 
bad  disadvantage  if  X2  is  expected  to  reach  zero  in  finite  time. 

A 

It  is  evident  that  the  steady  state  error  in  X2  or  x^  is  directly 
proportional  to  the  increase  or  decrease  in  parameters.  At  certain  values 
of  the  parameters  there  may  be  no  steady  state  error  at  all,  or  a small 
steady  state  error  and  for  large  deviations  in  parameters  a large  steady 
state  error.  The  compounded  effect  may  be  disadvantageous. 

In  a physical  system  if  the  parameter  deviations  are  small,  0.01 
percent,  and  if  small  steady  state  error  is  tolerated,  then  the  lower  order 
observer  will  serve  well.  If  not,  then  another  lower  order  observer  will 


tinve  to  be  designed.  Note  that  even  though  the  eigenvalues  selected  make 
the  estimated  states  decay  to  zero  in  2 seconds,  the  parameters  may  cause  a 
big  steady  state  error. 

Therefore  this  observer  may  or  most  probably  may  not  be  used  to 
estimate  the  two  states  X2  and  ■ 


3.3.  Example  No,  2 

From  the  equations  of  the  famous  stick-cart  system,  the  linearized 
stick  part  is  given  by 


Xi  = X2 


X2  = aX^  + bu 


where,  Xj^  is  linear  displacement  and  X2  velocity. 


L(m  +4M) 


L(m  +4M) 


wliere  M is  tlie  mass  of  the  cart 

m is  the  mass  of  the  stick 
1.  is  half  the  length  of  the  stick 
g is  the  acceleration  due  to  gravity, 
Typical  values  of  these  parameters  are: 

g°  » 9.78  m/sec^ 

M°  - 10  Kg 
m°  - 1.24  Kg 


(3.3.4) 


1 meter. 


'54 


wlicre 


Equation  (3.3.1)  can  be  put  into  the  form 


X = 


or 


o 


1 


0 

i1 

X + 

[ol 

3,, 

L 1 

2 

L . 

= _a  = - 3g°(m%M°)  ^ 

0,0  , ^ 

L (m  +4M) 


y = X, 


-8 


(3.3.5) 


(3.3.6) 


e 


o 

2 


= b 


3 

0,0  , o > 

L (m  +4M  ) 


-0.0727. 


Thus  the  eigenvalue  will  be  + v'^. 

Tn  this  system  there  are  four  parameters  and  they  are  given  as 

f o 1 1 aws . 


(3.3.7) 


All  tliose  parameters  are  subjected  to  slight  changes  and  thus 
could  affect  the  estimated  states. 

It  is  required  that  X2  be  estimated. 

The  Luenberger  observer  for  system  (3.3.1)  is  given  by 


Z = [-k]z  + L-k^+a°]y  (3.3.8) 

X2  - Z+ky.  (3.3.9) 

The  coefficient  k is  selected  to  be  more  negative  than  the  most  negative 


eigenvalue  of  (3.3.1). 


I 
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k V 8 


(3.3.10) 


H 

I 

I 


Using  ttic  same  techniques  as  adopted  in  Chapter  2,  we  design  the 
open-loop  sensitivity  model.  The  model  is  given  by  (3.3.11) 


0 


I -a" 

i_  1 


h-02 


Tlie  transformation  matrix  T and  are  given  as 


, = [o  i; 


(3.3.11) 


(3.3.12) 


I 

I 


As  there  are  four  parameters  and  one  estimated  state,  there  will  be 
four  sensitivity  functions.  TTiese  sensitivity  functions  are  given  below,  and 
they  can  be  generated  from  the  sensitivity  model  (3.3.11). 


= ( ) (m°+M°)h.  „ 


0.06h 


6X 


2 

6m 


I,  (m  +4M  ) 


O O 

27£_m_ 
o 


1,2 


m%4M°  2 


L“(m”+4M°)^ 

= -0.0047h^  2 ■ 0-097X2 


6X.  o o 

— s 2 7b  m . _ 1 V 

,q2,  o 3 "l,2  o.,o  2 

(m  +4M  ) m +4M 


0.038h^  2 ■ 0-024X2 


^ _ 9g°(m°  +M°)  - i X 


- -0.6h^^2  - ^2 


(3.3.13) 


(3.3.14) 


(3.3.15) 


(3.3.16) 


rrrrir3Bf*sac??!??^ 


We  only  require  !i.  „ out  of  the  three  known  sensitivity  functions  to  generate 
I-  > ^ 

tlic  sensitivity  functions  for  the  four  parameters  in  tliis  problem. 

Figure  1.6  is  the  diagram  for  generating  the  estimated  states 
sensitivity  functions  for  Example  2.  The  plots  of  the  sensitivity  functions 
(3.3.13)  to  (3.3.16)  are  illustrated  in  Figures  3.7  and  3.8. 

Common ts 

From  Figures  3.7  and  3.8  it  is  noticed  that  the  sensitivity  is 
in  the  order  of  10^^  and  therefore  the  observer  of  Example  2 is  highly 
sensitive  to  the  parameters  of  the  system.  This  is  a very  good  example  for 
the  type  of  systems  for  which  lower  order  observers  cannot  be  designed.  Thus 
the  observer  of  Example  2 is  useless,  as  it  will  make  tremendous  errors  in 
the  estimation. 


3.4.  Example  No.  3 

The  equations  of  a mercury  thermometer  from  [6]  are  given  by  the 


equation 


where 


*r  A- 

m g m g m 


— e 

R R o 
m g 


C 

8 


9 

m 


0 

o 


is  the  resistance  of  the  glass 
is  the  capacitance  of  the  glass 
is  the  resistance  of  mercury 
is  the  capacitance  of  the  mercury 
is  the  temperature  of  the  mercury 
is  the  temperature  of  the  environment. 


(3.4.1) 


40 


Typical  values  of  the  parameters  for  10  lbs.  of  glass  and  100  lbs. 

of  mercury 

= 1/5  Degree/(Btu/minute) 

= 2 Btu/degree 

® (3.4.2) 

r|^  = 1/60  Degree/ (Btu/minute) 

= 3 Btu/degree. 
m 

Equation  (3.4.1)  can  bo  put  into  the  form 


0 

“1 

1 

0 

-a 

-Of 

X + 

e 

1 

2 

2 

y = X, 


wliere 


o o„o,  o 
C R C R 
m m g g 


= 50 


= ( 


1 


o o 
C R 
m m 


o o 
C R 
g g 


o o 
C R 
g m 


•)  = 52.5 


1 


o o o o 
C R C R 
m m g g 


- 50. 


There  are  four  parameters  in  this  system.  They  are 


m 


®2  “ ^m 
«3  - ‘'g 
«4  ■ 


(3.4.3) 


(3.4.4) 


(3.4.5) 


(3.4.6) 


(3.4.7) 


41 


All  these  parameters  are  affected  by  the  conditions  of  temperature  and 
physical  materials,  tlierefore  they  tend  to  change.  A detailed  study  on  the 
topic  of  heat  transfer  will  verify  this  assumption. 

It  is  required  tliat  state  X2  be  estimated.  Note  that  X2  represents 
the  temperature  at  the  inner  surface  between  the  glass  case  and  mercury.  The 
l.uenbergcr  observer  for  system  (3.4.1)  is  given  as 


Z = [-(k+o'°)]z  + L -(k^ +a°)]y  + I'e^u 


X2  = Z + ky 


+k)  must  be  more  negative  than  -50. 


The  sensitivity  model  for  the  open-loop  sensitivity  functions  is 


given  by 


h-0"  X, 

^ ■ 


The  transformation  matrices  are  given  as 


T = P., 


1 0 


0 1 


, N « [O  l]. 


As  there  are  four  parameters  and  one  estimated  state,  there  will  be 
four  sensitivity  functions  which  can  be  written  as 

‘*2  1 - 

sr  ■ ■ - S.2  - ^ '‘2 

m m 


-833. 33h^  2 " 0-33X2 


A2 


Cv 

^ ^ h - X 

m m 

= -15  X lo\  - 60X 

I fZ  1. 


^ ^ _L' 

“ ■ C°  ''*2  - "2 

g g 

= -1250hj^  2 ■ 0-5X2 


ftX-  P°  , 

— — = - — h - — X 
5R  „o  "i,2  „o  ^2 

g R„  ’ R 
g g 

= -12*5  X lo\^  2 • 5X2- 


(3.4.13) 


(3.4.14) 


(3.4.15) 


These  are  the  four  sensitivity  functions  that  are  generated  from 
the  sensitivity  model  (3.4.10),  Figure  3.9  is  the  diagram  for  generating 
the  estimated  states  sensitivity  functions  for  Example  3.  The  plots  of  these 
sensitivity  functions,  (3.4.12)  to  (3.4.15)  are  illustrated  in  Figures  3.10 
and  3.11. 

Comments 

From  Figures  3.10  and  3.11  we  notice  that  all  the  sensitivity 
6X2 

functions,  except  » decay  to  zero  by  three  seconds.  The  parameter  R^, 

m 

or  the  resistance  of  mercury  causes  a steady  state  error  to  exist  in  the 

A 

temperature  X2  or  the  temperature  at  the  inner  surface  between  the  glass  case 
and  mercury.  A ten  percent  increase  in  the  resistance  will  cause  a 1.66 X 10 
deviation  in  the  temperature  X2'  temperature  itself  reaches  a steady 

-3 

state  of  -30 X 10  the  positive  deviation  is  not  considerable. 
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It  is  noticed  that  tlie  deviations  of  tlie  capacitances  C and  C are 

m g 

-3 

of  the  magnitude  10 

-3 

The  deviation  due  to  R is  also  of  the  magnitude  10  . The  maximum 

g 

or  minimum  deviations  are  negligible  as  compared  to  the  maximum  and  minimum 
values  of  tile  temperature  . It  must  he  mentioned  though  that  the  resistance 
of  mercury  causes  a slight  steady  state  error  or  deviation  in  steady  state. 

Thus  it  is  seen  that  the  lower  order  observer  for  Example  3 is 
reasonably  independent  of  variations  in  parameters  and  thus  is  a good  open- 
loop  estimator. 


3.5.  Example  No.  4 

A system  is  described  by  the  equations 


Xi  = a^X^ 

^2  ~ ^2^1  ^3^2 


(3.5.1) 


y = X, 


(3.5.2) 


where  a^^,  a^,  a^,  and  b are  assumed  to  be  parameters  that  change  very 
slightly.  The  nominal  value  of  these  parameters  are 


2 I 0^2^ 


-2,  a"  - -3,  b^ 


(3.5.3) 


There  are  four  parameters  and  they  are  given  by 


84  ■ *’• 


(3.5.4) 
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It  is  required  to  estimate  the  state  The  Luenberger  observer  for 

system  (3.5.1)  is  given  by 


Z = - kaj^  JZ  + I -(k  - a^k  - a2>  ly  + I b]u 


= X + ky. 


(3.5.5) 

(3.5.6) 


The 


eigenvalues  of  (a^-kaj^)  must  be  more  negative  than  the  eigenvalues  of 


3 v7 

(3.5.1),  that  are  - j — i — ^ • 


The  sensitivity  model  for  the  open-loop  sensitivity  functions  is 


given  by 


h“  = 


0 o 
-or,  -O' 

1 2 


ba^  1 


(3.5.7) 


where 


o 

"i  = 


° ° - / 
■^1^2  " 


-3. 


The  transformation  matrices  are  given  as 


(3.5.8) 


T = b 


, N = [O  iJ 


(3.5.9) 


The  sensitivity  functions  are  now  generated  from  the  sensitivity 
model  (3.5.7)  and  are  given  by 


6a, 


-ba°  h^^2  ■ ^ S,2 


(3.5.10) 


6a., 


-ba°  h^^2  * ^.2 


(3.5.11) 


4 a 


The  plots 
and  3.14.  Figure  3. 
for  E.xample  4, 
Comments 


'T  " '*2  ^ 


-2  !i 


2,2 


(3.5.12) 


5X 


ST  ■ f h ■ h' 


(3.5.13) 


of  these  sensitivity  functions  are  given  in  Figures  3.13 
12  is  tlie  diagram  for  generating  the  sensitivity  functions 


Of  the  four  examples,  this  example  has  the  best  lower  order  observer, 

as  far  as  steady  state  error  goes.  The  deviations  will  be  of  the  order  of 
-2 

10  . This  can  be  seen  in  Figures  3.13  and  3.14.  The  sensitivity  functions 

reach  steady  state  of  zero  at  an  average  of  6 sec.  The  influence  that  the 
different  parameters  have  on  the  estimated  state  can  be  easily  seen  in  the 

figure.  The  four  parameters  have  a different  effect  on  for  example,  a^^  has 
an  opposite  effect  as  does  a^ . They  both  cancel  each  others  effect  on 
Tlius  tlie  parameters  a^  and  b have  a greater  influence,  and  this  is  visible  in 
Figure  3.14.  Both  parameters  a^  and  b have  an  effect  on  the  overshoot, 
undershoot  and  the  steady  state. 

Thus  it  can  be  concluded  that  the  lower  order  observer  of  E.xample  4 
is  quite  independent  of  parameter  variations  and  can  thus  be  safely  called  a 


good  observer  or  an  estimator. 


T I I I I « I « I I I 
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Time,  in  sec 


I Figure  3.13  Kstimated  state  ^2  sensitivity  functions 

for  Example  U. 
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3.6.  Analysis  o£  the  Sunsitivity  of  Lower  Order  Observers  in  Feedback  Loops 
In  tlic  introduction  to  this  chapter  we  have  outlined  the  various 
points  to  be  investigated.  Some  of  these  points  fall  into  the  feedback  category 
and  will  be  analyzed  here  in  detail.  In  the  case  of  observers  in  feedback  we 
will  study  three  main  points  of  interest.  They  are;  general  sensitivity  or 
effect  of  variation  of  parameters  on  closed-loop  estimated  states,  effect  of 
eigenvalues  of  an  observer  on  sensitivity  and  finally  comparison  between  the 
lower  order  observer  and  the  full  order  observer  in  sensitivity.  To  accomplish 
this  task  an  experimental  analysis  will  be  undertaken.  Examples  5,  6,  and  7 
will  be  analyzed  at  the  end  of  tliis  section. 


3.6.1.  Example  No.  5 

A three-parameter  linear  time-invariant  system  is  given  by 


0 

1 

0 

1 

0 

! 

X = 

0 

0 

1 

X + 

0 

0 



-0 

2 

■‘^3 

I 

TP 

1 

(3.6.12) 


The  nominal  values  of  these  parameters  are  ^3  ~ ^ 

feedback  the  eigenvalues  must  be  placed  at  (-1,-3, -4). 

A lower  order  observer  and  a full  order  observer  is  required  to  be 
designed.  The  states  X2  and  have  to  be  estimated.  The  plots  of  the 
deviation  in  the  estimated  states  to  a ten  percent  increase  in  parameter 
value  must  be  made. 

In  order  to  determine  the  effects  of  change  of  eigenvalues  on  the 
sensitivity  of  the  lower  order  observer,  at  least  three  different  eigenvalues 


.1  CP—’ 


must  be  selected. 
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( 

f 


I 

I 

I 

\ 

? 

» 

I 

I 

I 
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Results : The  plots  of  these  deviations  for  eigenvalues,  (-8,-10), 

(-11,-13),  and  (-15,-20)  are  in  figure  3.15  to  Figure  3.20.  The  full  order 
observers  estimated  states  deviations  are  plotted  in  Figures  3.21  and  3.22. 

A A 

The  ma.ximum  and  minimum  deviations  in  the  estimated  states  and 
to  a ten  percent  increase  in  the  parameters  are  tabulated  in  Table  3.1. 


3.6.2.  Example  No.  6 

A tivo  parameter  system  is  given  by  the  equations 

u,  y = X^  (3.6.13) 

where  the  nominal  values  of  the  parameters  are 

= -8  and  0°  = 2 , 

It  is  required  to  estimate  X^ • The  feedback  must  be  designed  to 
adjust  the  eigenvalues  of  the  system  at  (-3,-4).  A similar  procedure  as 
adopted  in  Example  5,  is  followed  and  the  following  results  are  achieved. 

A 

Results : The  plots  for  the  deviations  in  X2  were  recorded  for  three  sets  of 

observer  eigenvalues  which  were  (-10),  (-20),  and  (-40).  The  plots  for  the 
full  order  observer  was  also  recorded.  These  plots  are  available  from 
Figure  3.23  to  Figure  3.26.  The  maximum  and  minimum  deviations  to  a ten 
percent  increase  in  parameters  were  tabulated  in  Table  3.2. 
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0 

X = 

0 

X 

+ 

3. 

3.6.3.  Example  No.  7 

A linear  time-invariant  system  is  given  by  the  equations 


0 

|0 


‘1 


X + 


X. 


(3.6.14) 
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Figure  3.15  Estimated  state  deviation  for  a ten  percent 
increase  in  paraneters , for  Example  5.  Observer 
eigenvalues  at  ( -8,  -10  ). 
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Figure  3.18  Estimated  state  deviation  for  a ten  percent 
increase  in  parameters,  for  Example  5.  Observer 
eigenvalues  at  ( -11,  -13  ). 
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Increase  In  parameters,  for  Example  5.  Observer 
eigenvalues  at  (-15,-20  ). 
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Figure  3.21  Estimated  state  X2  deviation  for  a ten  percent 

increase  in  parametes,  for  a Pull  Order  Observer 
of  fexample  5. 
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Tabic  3 . 1 
For  Kxamplc  5 

Max.  & Min.  Increment  for  a 10?  Increase  in 


I.ower  Order 

r 

Lower  Order 

Lower  Order 

Observer 

Observer 

Observer 

Full  Order 

Fi gen values  at 

Eigenvalues  at 

Eigenvalues  at 

Eigenvalues  at 

(-8,-10) 

(-11,-13) 

(-15,-20) 

(-8, -8, -10) 

- 

0.6  X 10'^ 

0.6  X lO"^ 

0.2  X lo’^ 

1.6  X 10'^ 

^X3  (0-2 ) 

L 

-1. 1 X 10'^ 

-0.6  X lo'^ 

-0.5  X lO’^ 

-1.3  X lo'^ 

p 

0.8  X 10'^ 

0.75  X 10"^ 

0.4  X 10'^ 

2.0  X 10'^ 

«,(»3) 

-0.85  X 10'^ 

-0.46  X 10‘^ 

-0.36  X 10'^ 

-0.9  X lO'^ 

p 

0.7X  10'^ 

0.6  X lo"^ 

0.1 X lo'^ 

2.0  X 10'^ 

5X3(33) 

— 

-1.2  X 10‘^ 

-1.39  X 10'^ 

-1.87  X 10'^ 

-0.6  X 10'^ 

p 

1.3  X 10'^ 

0.65  X 10'^ 

0.56  X lo'^ 

0.4  X 10'^ 

'X2  (0^2  ) 

_ 

-2.7  X lo'^ 

-2.6X  10'^ 

-1.48  X lo'^ 

-0.8  X lO'^ 

p 

0.4  X 10'^ 

0,40 X lo’^ 

0. 18  X lO"^ 

1.6  X 10'^ 

5X2  (‘^3) 

_ 

-0.8  X lO"^ 

-0.5  X lO"^ 

-0.38  X lO’^ 

-0.7x  10'^ 

1,61 X lo'^ 

1.8  X lo'^ 

1,68  X lO'^ 

0.42  X lo"^ 

3X3(33) 

— 

-0. 1 X lo”^ 

0.0 

0.0 

-1.2  X 10‘^ 

Time,  in  sec 

Figure  3.23  Estimated  state  X2  deviation  for  a ten  percent 
T increase  in  parameters,  for  Example  6.  Observer 

' eigenvalue  at  ( -10  ), 
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Figure  3.24  Estimated  state  deviation  for  a ten  percent 
^ increase  in  parameters,  for  Example  6,  Observer 

) eigenvalue  at  ( -20  ). 
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Figure  3.25  Estimated  state  deviation  for  a ten  percent 
Increase  in  parameters,  for  Example  6.  Observer 
eigenvalue  at  ( -40  ). 
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Figure  3.26  Estimated  state  X2  deviation  for  a ten  percent 


increase  in  parameters,  for  Example  6.  Full 
Order  Observer, 
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Table  3.2 
For  Example  6 

. & Min.  Increment  6X  for  a 10%  Increase  in  Parameters 


Lower  Order 
Observer 
Eigenvalues  at 
(-10) 

Lower  Order 
Observer 
Eigenvalues  at 
(-20) 

Lower  Order 
Observer 
Eigenvalues  at 
(-40) 

Full  Order 
Eigenvalues  at 
(-10,-15) 

0.0 

0.0 

0.0 

0.8  X 10“^ 

-2 

-2 

-2 

-\.0>  10  ^ 

-0.95  X 10 

-0.9  X 10 

-2.8  X 10 

-2 

-2 

-2 

, 1 

1. 1 X 10 

1.05  X 10 

1.0  X 10 

0.2  X 10 

0.0 

0.0 

0.0 

-0.84  X lo'*' 

fi8 


The  nominal  value  of  those  parameters  are 


o - o o T a 
= 2,  a^  = 3,  P 


2. 


aj  is  considered  to  be  a constant,  therefore  the  system  has  two  parameters, 
a^  and 

The  system  has  to  be  stabilized  by  state  feedback,  in  sucli  a way 
tliat  tlic  eigenvalues  of  the  system  are  set  at  (-3,-4).  The  state  vector 
has  to  be  estimated. 

The  procedure  as  adopted  in  Example  5 is  followed  in  this  example. 
Results : The  plots  of  the  deviations  of  the  estimated  state,  for  a ten  percent 
increase  in  parameter  values  is  recorded  in  Figures  3.27  to  3.30.  The  maximum 
and  minimum  deviations  for  three  sets  of  observer  eigenvalues,  (-17),  (-23), 
and  (-37)  were  tabulated  in  Table  3.3.  Also  the  maximum  and  minimum  deviations 
of  for  the  full  order  observer  were  tabulated  in  Table  3.3. 


3.6.4.  Effects  on  Sensitivity  by  the  Eigenvalues  of  the  Lower 
Order  Observer 

The  eigenvalues  of  the  lower  order  observer  have  an  effect  on  the 
sensitivity  of  the  lower  order  observer.  This  conclusion  is  drawn  from  the 
observations  made  in  Tables  3.1,  3.2,  and  3.3.  In  all  the  three  examples, 
it  has  been  noticed  that  there  exists  a pattern  in  the  deviations  due  to 
increasing  (negatively)  eigenvalues.  The  eigenvalues  considered  were 
allowable,  that  is,  they  fulfill  the  necessary  conditions  as  defined  for  the 


Luenberger  observer. 


Time,  in  sec 

Figure  3.27  Estimated  state  deviation  for  a ten  percent 
increase  in  parameters,  for  Example  7.  Observer 
eigenvalue  at  ( -23  ). 
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I 


increase  in  parameters,  for  Example  7.  Observer 
eigenvalues  at  ( -23  ). 
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Figure  3.29  Estimated  state  deviation  for  a ten  percent 
increase  in  parameters,  for  Example  7.  Observer 
eigenvalue  at  ( -37  ). 
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Figure  3.30  Estimated  state  X2  deviation  for  a ten  percent 
increase  in  parameters,  for  the  Full  Order 
Observer  of  Example  7. 
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Table  3.3 
For  Example  7 

Max.  & Min.  Increment  for  a 10’4  Increase  in  Parameters 


t>X2(32) 


Lower  Order 

Lower  Order 

Lower  Order 

Observer 

Observer 

Observer 

Full  Order 

Eigenvalues  at 

Eigenvalues  at 

Eigenvalues  at 

Observer  at 

(-17) 

(-23) 

(-37) 

(-15,-17) 

0.32  X lO'^ 

0.3  X 10‘^ 

0.28  X lo'^ 

2.0  X lo'^ 

-0.95  X lO'^ 

CM 

1 

O 

X 

o 

-0.85  X 10‘^ 

-6. OX  lO'^ 

1.5  X lO'^ 

1 .4  X lo'^ 

1.3  X lo"^ 

CSJ 

o 

V, 

CN 

-0.7  X lo"^ 

-0.65  X lO'^ 

-0.6x  lO'^ 

-1.35  X lo'*- 
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It  has  been  noticed  tliat,  as  the  eigenvalues  of  the  lower  order 
observer  become  more  negative  than  the  previous  eigenvalues,  the  deviations, 
nuiximum  and  minimum,  in  the  estimated  states  are  reduced.  Thus  it  is 
concluded  tliat  if  given  a set  of  allowable  eigenvalues,  the  most  negative 
of  tliese  eigenvalues  will  make  the  lower  order  observer  less  dependent  on 
the  parameters  of  tlie  system. 

3.6.5.  Comparison  Between  the  Sensitivity  of  the  Lower  Order 
Observers  and  the  Full  Order  Observers 
In  all  the  three  examples,  it  has  been  noticed  that  the  full  order 
observer  is  highly  sensitive  as  compared  to  the  lower  order  observer.  The 
Tables  3.1,  3.2,  and  3.3  are  a proof  to  this  effect.  The  magnitude  between 
the  sensitivity  of  the  two  types  of  observer  is  also  large.  Thus,  compelling 
us  to  conclude  that  the  lower  order  observer  is  less  sensitive  to  parameter 
changes  than  a full  order  observer.  Note  that  the  eigenvalues  of  both  the 
observers  were  designed  to  be  identical  or  very  close  to  the  other. 

This  concluded  our  experimental  results.  The  DEC-10  computer  was 
used  to  do  tlie  computation  needed.  The  type  of  problems  used  are  some  typical 
types  and  were  chosen  to  Illustrate  some  aspects  of  sensitivity  of  the  lower 


order  observers. 
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4 . CONCLUS IONS 

This  study  is  divided  into  two  parts,  Tlic  first  was  to  develop 
the  sensitivity  model  for  the  lower  order  observer,  and  that  was  carried 
out  in  Chapter  2.  The  second  part  was  to  analyze  the  sensitivity  of  the 
lower  order  observer.  The  experimentation  and  results  were  carried  out  in 
Chapter  3,  Now  we  shall  comment  on  the  total  study  and  on  the  subject  of 
sensitivity  of  lower  order  observers,  in  general.  It  must  be  noted  at  this 
point  that  our  conclusions  are  drawn,  in  the  case  of  the  sensitivity 
analysis  from  the  results  of  seven  examples  that  were  analyzed  theoretically 
and  with  the  help  of  the  DEC-10  computer. 

The  sensitivity  model  developed  is  a lower  order,  n-dimensional 
model.  In  certain  cases  the  dimensions  of  the  model  is  (n-r).  All  the 
sensitivitv  functions  are  obtained  from  this  model.  Thus  we  have  succeeded 
in  reducing  the  dimensions  of  the  model. 

The  sensitivity  of  the  lower  order  observers  has  given  some  insight 
into  the  design  of  these  observers.  The  observed  properties  are  the  following. 
Not  all  systems  can  have  a reasonable,  parameter  dependent  lower  order 
observer.  There  are  other  systems  that  can  have  lower  order  observers  only 
under  certain  conditions.  In  general,  systems  may  be  able  to  have  a reason- 
ably good  parameter  dependent  lower  order  observer  that  will  estimate  the 
states  of  the  system  without  excessive  errors  in  the  estimation.  Moreover, 
the  eigenvalues  of  the  lower  order  observer  of  a system,  effects  the  sensi- 
tivity of  the  lower  order  observers.  From  a set  of  eigenvalues  that  satisfy 
the  necessary  conditions  for  designing  the  parameter  independent  lower  order 
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observer,  the  most  negative  eigenvalues  will  make  the  lower  order  observer, 
least  dependent  on  the  parameters  of  the  system.  Therefore  when  design- 
ing a lower  order  observer  it  will  be  useful  to  keep  this  point  in  mind 
in  order  to  design  a good  estimator. 

The  lower  order  observer  is  less  dependent  on  the  system  parameters 
than  the  full  order  observers.  Therefore,  the  lower  order  observer  must  be 
used  wherever  possible  because  it  is  less  dependent  on  the  parameters  of  the 
system  and  therefore  will  estimate  better. 

This  concludes  our  study  of  the  sensitivity  of  the  lower  order 
observer.  All  the  conclusions  reached  are  due  to  some  theoretical  and 
some  experimental  results.  The  Luonberger  Observer  was  used  to  estimate 
the  states. 

It  must  be  mentioned  here  that  it  is  possible  that  the  experimentation 
on  other  models  of  lower  order  observns  may  deviate  from  the  conclusions 
reached  in  this  thesis.  There  is  no  proof  that  these  conclusions  are  true, 
but  experimental  results  have  shown  that  there  may  be  some  truth  to  our 
cone lus ions . 
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